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Abstract 

In  this  paper,  we  study  the  stability  and  evolution  of  the  solitary  waves  in  perturbed  generalized 
nonlinear  Schrodinger  equations.  Our  method  is  based  on  the  completeness  of  the  bounded  eigen¬ 
states  of  the  associated  linear  operator  in  L2  space  and  a  standard  multiple-scale  perturbation 
technique.  Unlike  the  adiabatic  perturbation  method,  ours  uncovers  all  the  instability,  mechanisms 
in  the  perturbed  equations.  As  an  example,  we  consider  the  perturbed  cubic-quintic  nonlinear 
Schrodinger  equation  in  detail  and  determine  the  stability  regions  of  its  solitary  waves.  The  gener¬ 
alization  of  this  method  to  other  perturbed  nonlinear  wave  systems  is  also  discussed. 
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1  Introduction 


In  recent  years,  the  perturbed  generalized  nonlinear  Schrodinger  equation  has  attracted  a  great 
deal  of  attention.  This  equation  is  of  the  form 

iAt  +  Axx  +  f{\A\2)A  =  ep(A,  Am),  (1-1) 

where  /  is  a  real-valued  algebraic  function,  p  is  a  spatial  differential  operator,  and  e  is  a  small 
parameter.  It  has  been  shown  to  govern  the  evolution  of  a  wave  packet  in  a  weakly  nonlinear  and 
dispersive  medium  and  has  thus  arisen  in  diverse  fields  such  as  water  waves,  plasma  and  nonlinear 
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optics  [1.  2.  3.  4.  5].  In  particular,  this  equation  is  now  widely  accepted  in  optics  field  as  a  good 
model  for  optical  pulse  propagation  in  nonlinear  fibers  (see  [6.  7|  and  the  references  therein).  The 
rapid  advances  in  optical-soliton  based  fast-rate  telecommunication  systems  in  recent  years  has 
stimulated  intensive  research  on  it.  Another  application  of  Eq.  (1.1)  is  in  pattern  formation,  where 
it  has  been  used  to  model  some  nonequilibrium  pattern  forming  systems  (see  (8|  and  the  references 
therein).  The  unperturbed  form  of  this  equation  (e  =  0)  supports  solitary  waves  of  the  form 

A  =  eiVx/2^'-v2/4)t-lpoa0(x  -  Vt  -  x0),  (1.2) 

where  a0{9)  is  a  real- valued  function  and  satisfies  the  equation 

flO 86  ~  +  f{a o)ao  =  0.  (1.3a) 

do  — *  0,  |(9|  — ►  oo,  (1.3b) 

and  V'. .j(>  0),  jo  and  po  are  arbitrary  real  constants.  Note  that  the  existence  of  these  solitary 
waves  is  the  basis  for  telecommunication  systems  using  optical  solitons  as  information  bits.  When 
perturbations  are  present,  one  very  important  concern  is  whether  these  solitary  waves  will  persist 
or  not.  This  question  has  been  studied  extensively  in  the  literature.  The  linear  stability  of  the 
solitary  waves  in  the  unperturbed  equation  (1.1)  has  been  investigated  in  (9,  10).  where  a  criteria  for 
instability  was  given.  When  f(x)  =  x,  Eq.  (1.1)  is  the  perturbed  nonlinear  Schrodinger  equation, 
whose  solitary  waves  were  examined  in  numerous  articles  such  as  [7,  11,  12,  13,  14,  15]  among 
others.  The  dynamical  behavior  of  the  solitary  waves  in  this  equation  is  now  well  understood. 
But  for  more  general  forms  of  /,  the  results  are  few  and  far  from  complete.  This  case  has  been 
investigated  in  [8,  16,  17,  18]  by  various  methods.  In  [8],  the  adiabatic  perturbation  technique  was 
employed  (see  also  (13]).  These  authors  assumed  a  quasi-stationary  form  for  the  solitary  wave, 
determined  the  slow  evolution  of  the  parameters  of  this  wave  and  then  discovered  certain  types  of 
instability  from  those  evolution  equations.  But  as  they  pointed  out,  the  stability  they  established 
only  refers  to  the  particular  class  of  perturbations  compatible  with  the  quasi-stationary  solution 
assumption.  It  was  recognized  that  there  could  be  other  instability  mechanisms  which  could  not  be 
found  by  this  adiabatic  method  (this  is  indeed  the  case).  In  (16],  the  Evans  function  approach  was 
used.  By  calculating  the  small  eigenvalues  bifurcating  from  the  zero  eigenvalue  of  the  associated 
linear  operator  in  the  unperturbed  equation  (1.1),  the  author  gave  the  conditions  for  this  type 
of  instability.  His  results  are  basically  equivalent  to  those  in  [8],  We  would  like  to  emphasize 
here  that  both  results  in  [8]  and  [16]  missed  certain  types  of  instability  hidden  in  Eq.  (1.1).  The 
numerical  approach  to  this  problem  was  taken  in  (17,  18].  In  this  work,  the  authors  investigated 
the  stability  of  analytic  solitary  waves  of  the  cubic-quintic  complex  Ginzburg- Landau  equation  and 
found  that  they  are  generally  unstable,  except  in  a  few  special  cases.  The  instability  was  caused  by 
the  existence  of  growing  disturbances  whose  largest  growth  rates  were  numerically  estimated.  The 
authors  also  obtained  parameter  regions  in  which  stable  solitary  waves  exist  for  various  choices  of 
parameter  values.  In  their  work,  e  was  not  small,  and  it  was  not  clear  just  what  was  the  source  of 
their  instability. 

In  this  paper,  we  develop  a  new  analytical  method  for  studying  the  stability  and  evolution  of  the 
solitary  waves  in  Eq.  (1.1).  This  method  is  based  on  knowing  the  closure  of  the  bounded  eigenstates 
of  the  associated  linear  operator  in  Li  space,  combined  with  a  standard  multiple-scale  perturbation 
method.  Li  is  the  space  of  all  the  square-integrable  functions.  In  essence,  this  method  is  similar 
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to  the  one  developed  in  [12]  (see  also  [14.  15))  for  solitons  in  perturbed  integrable  equations.  But 
here  the  new  feature  is  that,  since  the  unperturbed  equation  (1.1)  is  non-integrable  in  general,  the 
completeness  of  the  bounded  eigenstates  (or  equivalently,  the  Green  s  function)  of  the  associated 
linear  operator  has  to  be  established  anew.  We  will  use  a  direct  scattering  technique  analogous  to 
that  in  ;20.  21]  to  accomplish  this  task.  In  this  process,  the  structure  of  the  spectra  of  this  linear 
operator  will  also  be  obtained  and  detailed.  Using  this  new  method,  we  can  uncover  all  instabilities 
of  the  solitary  waves  in  Eq.  (1.1)  and  give  a  complete  account  of  the  stability  and  evolution  of  the 
solitary  waves  of  this  equation.  We  would  like  to  point  out  here  that,  in  principle,  this  method 
can  be  applied  to  any  perturbed  nonlinear  wave  equation  for  uncovering  all  the  instabilities  of  its 
permanent  waves.  We  will  come  back  to  this  point  in  section  4. 

After  the  general  procedure  of  this  method  is  introduced,  we  will  apply  it  to  the  perturbed  cubic  - 
quintic  nonlinear  Schrodinger  equation  and  carry  out  the  analysis  in  detail.  Assuming  the  pertur¬ 
bation  to  only  contain  terms  of  the  Ginzburg-Landau  type  (as  in  [8,  16.  17.  18]),  we  will  show  that 
the  perturbed  cubic-quintic  nonlinear  Schrodinger  equation  allows  at  most  two  solitary  waves,  of 
which,  at  most  one  is  stable.  We  also  find  that  the  solitary  waves  of  the  model  equation  (1.1)  have 
three  instability  mechanisms  which  are  related  to  perturbations  of  respectively  the  zero,  non-zero 
(discrete)  and  continuous  eigenvalues  of  the  associated  linear  operator  in  the  unperturbed  equa¬ 
tion  (1.1).  Of  these  instabilities,  the  instability  related  to  perturbations  of  the  non-zero  discrete 
eigenvalues  has  never  been  studied  before.  Its  capture  requires  expansion  of  the  perturbation  series 
of  the  solution  out  to  second  order,  e2.  We  further  derive  the  necessary  and  sufficient  stability 
conditions  for  these  solitary  waves  and  specify  the  regions  of  parameter  space,  inside  of  which,  this 
equation  has  stable  solitary  waves.  Finally,  the  generalization  of  this  method  to  the  study  of  the 
stability  of  permanent  waves  in  other  nonlinear  wave  systems  is  discussed. 


2  The  Procedure 


In  this  section,  we  detail  the  procedure  for  studying  the  stability  and  evolution  of  the  solitary  waves 
(1.2)  in  Eq.  (1.1).  For  simplicity,  we  consider  the  case  where  the  perturbation  term,  p,  is  of  the 

form  ”  SM 

p(a.a')  =  E  wd-4!2)  ajv  (2a) 

where  p*  (Ac  =  1, . . .  ,n)  are  complex  functions.  This  will  exclude  parametrically  forced  perturba¬ 
tions  (see  [22|).  But  even  in  such  cases,  the  analysis  given  here  can  be  readily  modified.  Anticipating 
the  slow  evolution  of  the  free  parameters  of  the  solitary  wave  when  a  perturbation  is  present,  we 
write  the  solution  of  this  equation  in  the  form 

A  =  eiV6/2+lf,a(e,t,TuT2 e),  (2-2) 


where 


ft  ft  v2 

9  =  x-  Vdt-Oo,  P=  (w  +  — —  )dt  -  po, 
Jo  Jo  4 


(2.3) 
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and  V.Oq  and  pn  are  all  functions  of  slow  time  Ti  -  tt. 
substituted  into  Eq.  (1.1).  the  equation  for  a  is  found  to  be 


f2t .  When  Eq. 


Kit 


jja  —  aoe  —  f(\a\2)a  — 


f.F  -  -  i0QTxae 

—f2{lQT2  —  i0oT2aQ  ^ 


\  y,  o  ,  ^ 

- f  T  Pr,Tx  )af 

^  -A)72)a}  -OO.3). 


i  2.2:  i? 


(2.4: 


where 

F  =  p(A..4*)e-‘19  2_tp. 

To  solve  this  equation,  we  expand  a  into  a  perturbation  series 

a  =  ao(0)  +  eai  +  e2Q2  +  •  •  • 


C2.5) 


(2.6) 


and  take  a0  to  satisfy  Eq.  (1.3).  Thus,  the  zeroth  order  of  (2.4)  is  now  trivially  satisfied.  At  order 
f f,  a\  is  governed  by  the  linear  equation 


iau  —  -va !  +  am  +  r{0)ai  +  q{9)a\  =  u*i . 
a  1 1 1  =o  =  0. 


(2.7a) 

(2.7b) 


where 


r  =  /(a2)  4"  ao /  (ao)’  'I  —  ao/  (ao)’ 


w\  =  Fq  —  iaoTi  +  Mori  aoe  ( 


ve. 


0T\ 


vT[e 


+  POTi  )<*0. 


(2.8) 

(2.9) 

(2.10) 


2  2 

F0=p(A0^)e-lV9/2-ip, 

and  40  =  eive/2*ipa0.  Note  that  F0  appears  to  have  a  fast  time  t  dependence,  but  it  actually 
does  not.  due  to  the  form  of  p  in  Eq.  (2.1).  This  fact  wiU  be  used  in  the  later  analysis.  Denoting 
Ux  =  ( a(.a;)r ,  Eq.  (2.7)  can  be  rewritten  as 


(idt  +  L)U{  =  (w\,  —iv[)T, 

CAIt-o  =0. 


where  the  linear  operator  L  is 


L  —  <73 


dee  ~  ^  +  r  q 

q  dee  ~  •*>  +  r 


(2.11a) 

(2.11b) 

(2.12) 


and 


o  i  = 


<73  = 


(2.13) 


are  Pauli  spin  matrices. 

The  linear  operator  L  is  the  key  to  solving  this  problem.  If  we  could  expand  all  quantities  in  the 
eigenstates  of  this  operator,  then  we  can  expand  and  solve  (2.4)  to  all  orders.  So,  first  we  need  t  e 
eigenstates.  w(x.  A),  where  A  is  the  eigenvalue,  and  the  eigenvalue  spectra  of  the  operator  L.  where 

Lip  -  A  ip.  (2-14) 
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We  will  now  discuss  the  main  features  of  this  eigenvalue  problem.  Due  to  the  form  of  L.  if  A  is 
an  eigenvalue  and  v  the  corresponding  eigenfunction,  then  it  follows  that  -A.  A'  and  -A'  are  also 
eigenvalues,  with  corresponding  eigenfunctions  as  (T\v.  l'*  and  cric  .  Ln  the  appendix,  we  detail 
the  structure  of  the  spectra  of  the  operator  L.  using  the  direct  scattering  technique.  We  show  that 
the  discrete  eigenvalues  of  L  are  zeros  of  the  analytical  function  A2  defined  there.  Now.  there  are 
four  free  parameters  in  the  unperturbed  solitary  wave  (1.2).  Perturbations  of.  or  shifts  in  these 
four  free  parameters  correspond  to  the  four  degrees  of  freedom  represented  by  the  degenerate  A  =  0 
eigenvalue.  Thus  A  =  0  is  at  least  a  four-fold  eigenvalue  of  L.  Two  of  these  degenerate  eigenvalues 
have  two  discrete  eigenfunctions 

t’oi  =  aoe(h  I)7”-  ©02  =  oo(T  —  l)r.  (2.15) 

which  exactly  satisfy  the  eigenvalue  equations 

Lrv oi=0,  L1V02  —  0-  (2-16) 

The  other  two  eigenvalues  correspond  to  two  generalized  eigenfunctions,  called  “derivative  states" 
in  [14.  15|  and  elsewhere,  and  are  given  by 

©01  =  2^ao(1-  -1)r’  ©02  =  aou)(T  1)T-  (2-17) 

These  states  are  not  true  eigenfunctions,  but  they  are  necessary  for  closure  (23).  They  satisfy  the 
modified  eigenvalue  equations 

Lc&oi  =  t/'oit  fi</> 02  =  ©02-  (2.18) 

The  total  number  of  the  discrete  eigenvalues  of  L  (with  the  multiplicity  of  all  degenerate  eigenvalues 
included)  is  given  by  the  “winding  number”  of  A2,  as  its  argument  moves  along  the  path.  P. 
described  in  the  appendix,  and  shown  in  Fig.  3.  The  continuous  eigenvalues  of  L  are  found  along 
the  two  half  lines  {A  :  A  >  <a>}  and  {A  :  A  <  -u>}.  For  each  value  of  A  in  those  intervals,  there  are  two 
continuous  eigenstates,  one  symmetric  in  9  (denoted  as  ws(9,  A))  and  the  other  one  anti-symmetric 
in  9  (denoted  as  t/>a(0,A)).  In  the  appendix,  we  also  show  that  the  union  of  the  discrete  regular 
eigenstates,  the  discrete  generalized  eigenstates  and  the  continuous  eigenstates,  form  a  complete 
set  in  the  space  of  L2.  Therefore,  Eq.  (2.11)  can  be  solved  by  expanding  the  solution.  U\ ,  and  the 
inhomogeneous  term,  in  this  complete  set  of  functions. 

For  simplicity,  in  the  rest  of  the  paper,  we  assume  that  zero  is  a  four-fold  discrete  eigenvalue  of  L. 
In  addition,  we  assume  that  L  has  only  two  other  simple,  discrete,  non-zero  eigenvalues,  denoted 
as  Ad  and  -A^,  with  the  corresponding  eigenfunctions  denoted  by  rbd  and  t-d  =  & i©d-  This  is  the 
case  for  the  perturbed  cubic-quintic  nonlinear  Schrodinger  equation,  to  be  discussed  in  more  detail 
later  in  this  paper.  If  L  has  more  than  the  above  eigenvalues,  or  if  the  non-zero  discrete  eigenvalues 
are  not  simple,  then  the  following  analysis  can  be  easily  appropriately  modified. 

Under  the  above  assumptions,  we  then  can  expand  (wi,-w[)T  and  U\  in  this  closed  set.  Using 
arbitrary  coefficients,  we  can  take 

(uq,  —w[)T  =  coii/>oi(0)  +  ^02^02(0)  +  doi0oi(0)  +  do2<£o2(0)  (2.19) 

+CdTpd{0)  +  C-d^-di^)  +  //{ca(^)1(,a($!  +  Cs(A)'0s(0,  A)}dA, 
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U\  =  hoivoi{9)  +  ho2vo2{0)  +  goiOoi{9)  +  502002(0) 

+hdUd{Q)  +  h-dV-d(Q)  f[{ha(^)va{9-  A)  +  hs(\)us(Q'  A)}dA. 

where  the  interval  l  =  (-oc.  -~’i  U  U.  oc).  We  will  define  an  inner  product  by 


(2.20) 


<f.g>=J  fT{e)<T3g(0)de.  |'2.21) 

Then  it  is  easy  to  show  that  the  only  non-zero  inner  products  of  these  bounded  eigenstates  are 
<  C'Ol.Oot  >•  <  t’02>  002  <  fi'd:  L'd  >■  <  C'-d,  t'-d  >.  <  I'a-^'a  >■  an<l  <  L's  >  • 


In  particular. 

<  b’a(--  A).  c'a(-.  A')  >=  ka(X)6(\  —  A'). 

<  t’3(-,  A).  v3(-.  A')  >=  MA)e(A-A'). 

where  jfca(A)  and  ks{\)  can  be  related  to  the  scattering  data  of  Eq.  (2.14).  When  the  expansions 
(2.19)  and  (2.20)  are  substituted  into  Eq.  (2.11).  and  the  above  inner  products  are  used,  then  we 
obtain  the  following  equations  for  the  coefficients  in  0\: 


(2.22a) 

(2.22b) 


.8  hoi  dho 2 

j— — f  501  =  On.  +  0°2  —  c°2’ 

.3501  .  ,%>2  , 

'-dT  =  d°'’ 


■dhd  ,  .  _ 

2  -1-  Qii 

at 


dh-d  x  l 

t  ^  A dh— d  C-di 

dh3 


i— - - b  A/ia  —  ca,  b  A/is  —  c„ 

ot  at 

hoi  =  ^02  =  501  =  502  —  tid  —  h-d  =  ha  =  h3  —  0.  at  t  =  0. 


(2.23a) 

(2.23b) 

(2.23c) 

(2.23d) 

(2.23e) 


C01 

< 

(lOl,  — u;i)r:Ooi  > 

C02 

<  (uq,  -w')T, 

CM 

0 

0 

> 

<  001,001  > 

<  002,  002 

> 

dot 

< 

(tui,  -u;J)r,0 01  > 

do2 

<  (wi, -w[)T , 

U’02 

> 

<  *001 ,  ^>01  > 

<  O’02  ,  002 

> 

Cd 

< 

(wi,  -w')T,Vd  > 

<  {w\.-w[)T.  1 

> 

<0d,0d> 

c-d  : 

<  0_d,  1 p-d 

> 

(2.24a) 

(2.24b) 

(2.24c) 


Ca  = 


<  (u’l,  -tt’*)r,P’a  > 

ka 


<  (u’l.  u.’*  )7~ .  dq  > 


Cs  = 


(2.24d) 


Note  that  c_d  =  —  since  0_d  =  <7id>d  and  is  real.  For  the  same  reason,  h_d  —  hj.  Since  u.'i  does 
not  depend  on  the  fast  time  (,  the  quantities  in  Eqs.  (2.24)  do  not  either.  To  suppress  the  secular 
terms  in  hoi)ho2i5oi  and  502 1  we  need  to  require  that 


coi  =  C02  =  dot  =  do2  —  0. 


(2.25) 
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Ln  view  of  Eqs.  (2.24)  and  (2.9).  these  four  conditions  will  produce  the  following  slow  evolution 


equations  for  V.jJ.9q  and  po  on  2~i  time  scale: 

*L  -  (2.26a) 

dTi  jrx  a20  dd 

± ±  =  2/j^  a0  Im(F0)rffl  .,gb) 

dTx  JTMhdB 

f"  a20d9  +  ^r  I"  =  2  /~  9a0  Im (F0)dO,  (2.26c) 

dl  \  7-oo  dl\  7 -x)  7-dc 

(»■£  +  2^)  r  (alLM  -  r  Halide  =  4  T  ^.Re(Fo)d9.  (2.264) 

aT\  dl  1  7-oc  Oil  7  -  00  7  -  00 


Here  "Re”  and  Tm”  represent  the  real  and  imaginary  parts  of  a  complex  number.  It  is  noted 
that  Eqs.  (2.26a.  b)  have  been  obtained  before  by  the  adiabatic  perturbation  method  in  18] . 
Similar  equations  were  also  derived  for  solitons  in  perturbed  nonlinear  Schrodinger  equations  (see 
(7.  11.  15]).  In  order  for  the  solitary  waves  of  the  model  (1.1)  to  be  stable,  these  equations  must 
have  stable  fixed  points.  Otherwise,  an  instability  will  arise.  Such  an  instability  would  be  due 
to  the  zero  eigenvalue,  of  the  linear  operator  L.  bifurcating  and  moving  into  the  unstable  region, 
because  of  the  perturbations.  This  has  been  discussed  in  [8,  16]. 


When  conditions  (2.26)  are  satisfied,  solving  Eq.  (2.23).  we  get 


hoi  —  ho2  —  3oi  =  302  =  0. 

(2.27a) 

hd  =  cd{l  -  ocd(T\)elXdt} / \d,  h_d  =  hmd. 

(2.27b) 

ha  =  ca{l-<*a(Tl)eiXt}/\, 

(2.27c) 

h3  =  cs{l  -  as(Ti)elAt}/ A, 

(2.27d) 

and 

(2.28) 

ad(0)  =  Qo(0)  =  as(0)  =  1. 

Then  the  solution  U\  is 

U\  =  hdWd(9)  +  h-dip-d{9)  +  j {ha{\)Va{9 ,  A)  +  hs(\)wa(9.  A)}dA. 

(2.29) 

Here  the  a's  are  constants  of  the  integration,  and  possibly  functions  of  7\,  as  indicated.  The  c  s 
are  slowly  varying  with  T\  when  V  and  u>  are. 


It  is  important  to  realize  here  that,  in  order  for  the  solitary  wave  (1.2)  to  be  stable,  in  addition  to 
the  conditions  (2.25),  we  also  need  to  require  that  the  coefficients  hd,h-d,ha  and  hs  in  do  not 
grow  unbounded,  on  either  the  t  or  the  T\  scales.  On  t  scale,  ha  and  h3  are  already  bounded  since 
the  continuous  eigenvalues  of  the  operator  L  are  always  real  (see  Eq.  (2.27  c,  d)).  But  in  order  for 
hd  and  h-d  to  remain  bounded  on  this  scale,  it  is  necessary  for  A^  to  be  real.  If  this  is  not  so,  then 
the  solitary  wave  is  unstable.  We  now  shall  assume  that  A<*  is  real. 
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On  the  T\  scale,  we  need  to  ensure  that  a oca  and  in  Eq.  (2.27)  remain  bounded.  To  obtain  the 
evolution  equations  for  these  coefficients,  we  need  to  expand  (2.4)  out  to  second  order,  a  .  When 
Eq.  (  2.6)  is  substituted  into  Eq.  (2.4)  and  terms  of  order  collected,  an  equation  for  <io  will  be 
obtained.  Denoting  U2  =  (<22.  ^2)^ •  equation  for  L‘2  is 


( idt  +  L ) 62  =  (u'2.  -w’2)T . 

(2.30a) 

6'2|t=o  =  0. 

(2.30b) 

where 

W2  =  F\  —  iaiTi  +  ^oTi aie  ~  (V90Tl/2  -  VTl9/2  +  P0T1  )a\ 

-iaoT2  +  i9oT2aoe  —  (W?or2/2  -  ^t29  / 2  +  Por2)ao 
-a0f,(al)al(al^2a[)-alr(al)(al+al)2/2. 

(2.31) 

F 1  =  {p.4(Ao,  .4^)Ai  +  p.4-(Ao.  .4^)A*}e  ,ve/2  ip , 

(2.32) 

and  Ai  =  elV6l2^lpa\.  This  equation  can  be  solved  analogous  to  Eq.  (2.11).  We  expand  (u;2,  -u'-2)r 

and  if2  as 

(u?2, 

—  U>2)T  =  Goi^oi(^)  +  0)2  V02  {9)  +  d^\0Q\{9)  +  d<)2<P02(9) 

+cdiPd(9)  +c-dw-d(9)  +  fl{ca{\)wa(9,  A)  +  cs(A)ips(0.  A)}dA. 

(2.33) 

and 

C/2  =  hoi  up  1  (9)  +  ^02^02  (9)  +  goi<Poi(9)  +  go2<Po2(9) 

Jrhdxi>d{9 )  +  h-dip-d{9)  +  j [{ha(\)xi>a{9 ,  A)  +  h3(\)U’s{9.  A)}dA. 

(2.34) 

The  coefficients  in  t/2  are  governed  by  equations  similar  to  (2.23)  with  only  a  hat  added  to  each 
quantity.  To  suppress  the  secular  terms  in  hoi,  ho2>  poi  and  302,  we  will  obtain  the  evolution 
equations  for  the  parameters  V,  jJ,  9q  and  po  on  T2  timescale.  These  equations  are  negligible  in 
deference  to  Eqs.  (2.26)  which  govern  their  evolutions  on  T\  timescale.  The  coefficient  hd  in  t/2  is 

governed  by  the  equation 

■  dhd  - 

^  0,  T  A d^d  Cdi 

at 

(2.35) 

where 

<  (u;2,  -ujj)7,^  > 

cd  — 

<  U-'d,  1 Pd  > 

(2.36) 

Now  dd  has  resonant  terms  which  are  proportional  to  elXdt.  To  see  this,  we  put  Eq.  (2.1)  into  (2.32) 

and  get 

Fl  =  +  ^loAT)  + 

(2.37) 

When  Eqs.  (2.27)  and  (2.29)  are  substituted  into  Eq.  (2.37),  we  find  that  the  e,A*£  and  e  ,Ai£ 
coefficients  in  F\  are  proportional  to  cdad/ and  (cdad)’  /  Xd  respectively.  Suppose  such  terms  in 

F\  are 

^ei(0)e<Ai£  +  (^e2(0)|*e-*A^£. 

A  d 

(2.38) 

then  the  coefficient  of  the  e,A<lt  term  in  is  found,  from  Eq.  (2.36),  to  be 

K  =  {i(cdad)ri  +  (k\  +  h2  +  kz)cdad}  / \d, 

(2.39) 
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where 


k[  =  (VdoTi /2  -  t  Pqt\) 


irjv2 


d  1 


'd2 


)d0 


A:  2  = 


-  vj2)de- 

(Cd  4-cj)  i'd2){2ao/'(ao)(t-'rfi  +  ^1^2  +  +  ao/"(ao)(L'rfi 

“IT"  I-Jvj  1  -  L&)d* 


(2.40a) 

^2)“ 


*3  = 


/T*.  ( &  1  t'd  1  +  e2t'd2  )d8 

LXM-^  ' 


(2.40b) 

(2.40c) 


and  ei  and  62  were  introduced  in  (2.38).  Since  KelXdt  term  in  c d  is  a  homogeneous  solution  of  Eq. 
(2.35),  then  in  order  to  suppress  the  secular  growth  in  hd<  we  must  have  K  =  0.  This  gives  us  a 
slow  evolution  equation  for  which  is 


£^l  =  l{kx+k2  +  k3)cdad.  (2-41) 

Note  that  both  k\  and  are  real  quantities.  Thus  if  \m(kz)  is  negative.  cdad  will  exponentially 
grow,  and  the  solitary  wave  (1.2)  will  be  unstable.  This  instability  is  caused  by  the  pumping  of 
the  solitary  wave  energy  into  the  discrete  eigenmodes  u'±d  when  the  perturbation  is  turned  on. 
It  can  also  be  interpreted  as  the  initially  real  discrete  eigenvalues.  moving  into  an  unstable 

region  in  the  presence  of  perturbations.  In  fact,  Ad  is  shifted  to  Ad  +  e(k\  +  £2  +  k%)  in  view  of  Eq. 
(2.27b).  This  instability  was  never  analyzed  before  in  the  literature.  Similar  argument  applies  to 
the  coefficients  ha  and  h3  of  the  continuous  eigenstates  in  C/2.  Suppression  of  the  secular  terms  in 
those  coefficients  will  produce  evolution  equations  for  aa(7\,  A)  and  as(Ti,  A)  on  T\  scale.  If  aa  or 
a3  grows  unbounded,  instability  will  also  arise.  This  instability  is  caused  by  energy  being  injected 
into  the  continuous  eigenmodes  of  L,  or  the  continuous  eigenvalues  of  L  moving  into  the  unstable 
region,  under  perturbations.  In  summary,  by  studying  the  evolution  equations  of  the  solitary  wave 
parameters  and  the  coefficients  in  the  U\  solution,  all  the  instability  mechanisms  can  be.  and  have 
been  now  uncovered. 


In  the  above,  we  have  obtained  the  slow  evolution  equations  for  V,uj.0q,po  and  ad.  The  equations 
for  aa  and  a3  in  the  coefficients  of  the  continuous  eigenmodes  of  U\  are  more  troublesome.  The 
reason  is  that  these  equations  involve  convolutions  which  couple  together  aa  and  a3  over  all  the 
continuous  eigenvalues.  To  circumvent  this  difficulty,  it  is  helpful  to  view  this  type  of  instability  as 
due  to  the  continuous  eigenvalues  of  L  moving  into  the  unstable  region.  Note  that  the  continuous 
eigenvalues  of  the  linearization  operator  around  a  solitary  wave,  even  for  the  perturbed  equation 
(1.1)  as  well  as  for  the  unperturbed  version,  can  easily  be  specified  (see  (24]) .  Thus  this  type 
of  instability  can  be  determined  without  the  necessity  of  deriving  and  examining  the  evolution 
equations  for  a<,  and  as.  With  this  hurdle  removed,  then  our  procedure,  as  described  above,  can 
be  carried  out,  and  the  stability  regions  of  the  solitary  waves  under  perturbations  can  be  specified. 
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3  The  Perturbed  Cubic-Quintic  Nonlinear  Schodinger  Equation 


Ln  this  section,  we  use  the  perturbed  cubic-quintic  nonlinear  Schodinger  equation  of  Ginzburg- 
Landau  type  as  an  example  and  carry  out  the  detailed  analysis.  This  equation  is  of  the  form 

iAt  +  Axx  -r  c3jA|2A  -r  e5iA|4A  =  ei( bxAxx  +  *  A  —  63|AjaA  -  65lA!4A).  C.iA) 


where  all  the  coefficients  are  real- valued,  and  c3  =  ±1  by  scaling.  When  e  =  0.  Eq.  (3.1)  supports 
solitary  waves  of  the  form  (1.2),  where 

|  i '  2  |  2  j 

a°(^  ^ c3  +  VT  +  16c5u;/3  cosh(2yC;  6) 

[f  c3  =  1.  this  wave  exists  when  c3u;  >  —3/16:  if  c3  =  —1.  it  exists  when  c3  >  0.  The  linear  operator 
L  is  given  by  Eqs.  (2.8)  and  (2.12)  where  f(x)  =  c3x  +  c5x2. 

We  first  establish  the  spectrum  structure  of  L.  We  have  shown  that  its  continuous  eigenvalues  are 
the  intervals  /  =  (-oc.  -jj\  U  [u>.  oo).  To  determine  the  total  number  of  its  discrete  eigenvalues,  we 
numerically  calculated  A2  along  the  path  P  shown  in  the  appendix  for  c3  =  ±1  and  c5,  w  being 
allowed  various  values.  For  each  of  the  three  cases:  (1)  c3  =  I.C5  >  0,  (2)  c3  =  l.c3  <  0.  and  (3) 
c3  =  -l  (c5  >  0),  the  results  are  always  qualitatively  the  same.  In  case  (2),  the  orbit  of  A2.  as  C 
moves  along  P.  is  sketched  in  Fig.  1.  We  see  that  in  this  case,  the  winding  number  of  A2  is  four, 
thus  L  has  four  discrete  eigenvalues  (multiplicity  of  eigenvalues  included).  In  the  other  two  cases, 
we  find  that  the  winding  number  of  A2  is  six.  Recall  that  A  =  0  is  always  a  discrete  eigenvalue  of 
L.  To  determine  its  multiplicity,  we  chose  a  small  closed  path  around  A  =  0  (i.e.  ^  =  e'*'4)  and  find 
that  the  winding  number  of  A2  is  always  four,  for  all  three  cases.  This  means  that  A  =  0  is  always 
a  four-fold  discrete  eigenvalue  of  L.  We  then  conclude  that  in  case  (2),  A  =  0  is  the  only  discrete 
eigenvalue  of  L.  while  in  the  other  two  cases,  L  has  two  additional  non-zero,  discrete  eigenvalues. 
Due  to  the  symmetry  of  the  eigenvalues,  these  two  non-zero  eigenvalues  have  to  be  either  real  or 
purely  imaginary.  In  addition,  one  is  always  the  negative  of  the  other.  We  will  denote  them  as  A d 
and  -Xd  as  before.  Closer  examination  reveals  that  in  case  (1),  A^  is  real,  and  in  case  (3),  it  is 
purely  imaginary.  This  is  consistent  with  the  results  in  [9,  10|.  It  indicates  that  in  case  (1).  the 
solitary  wave  (1.2)  is  linearly  neutrally  stable  in  the  unperturbed  equation  (3.1).  But  in  case  (3). 
it  is  linearly  unstable,  thus  also  unstable  under  weak  perturbations.  In  the  rest  of  the  paper,  we 
assume  c3  =  1. 


Next  we  determine  the  slow  evolution  equations  for  the  solitary  wave  parameters;  V,a;.#o  and  po ■ 
In  view  of  the  perturbation  term  in  Eq.  (3.1).  and  after  some  algebra,  we  find  that  Eqs.  (2.26) 


become 


dV_  =  _4b1f^a 
dTi  IZo  aoM  ’ 


(3.3) 


dw_  =  Sjyj)  —  b\V^J2  /3  4) 

&ln 

^2.  =  =  0  (3.5) 

dT\  dT\ 
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where 


3£>5-\  .  M  ,,  %5J'*x4d6 

S(*l>)  —  (27  —  b\jJ  )">■(.  2^3  t  ,  )  7x  2ja 

C5  4  4 c5  J_ x  a0da 

The  fixed  points  of  Eqs.  (3.3)  and  (3.4)  are  1=0  and 

S(ut)  =  0. 

[n  order  for  the  fixed  points  to  be  stable,  we  need  to  require  that 

s'M 


&1  >  o. 


£  In  f™xa'od0 


<  0. 


(3.6i 


(3. 


(3.8) 


Explicit  expressions  for  S(jj)  and  f^a^dO  can  be  obtained  depending  on  the  sign  of  c5. 


1.  c, 5  >  0:  In  this  case. 


where 


/  a^dO  =  y'3/cs  (tr/2  -  arctanu  l). 
y  ~  oc 

3  S3U 

~  16c?  51  +  S2UJ  +  tt/2  -  arctanu-1 


si  =  b\Cs  —  86305  +  965  —  327c|/3, 
s2  =  16c5(65  +  6^5/3), 


(3.9) 

(3.10) 

(3.11a) 

(3.11b) 


S3  =  —  61 C5  +  863  C5  —  965, 


(3.11c) 


and  _ _ 

u  =  y/l6c5u;/3.  (3.12) 

It  is  easy  to  check  that  S"(aj)  does  not  change  sign  for  uj  >  0.  hence  the  concavity  of  S(aj) 
does  not  change.  As  a  result.  Eq.  (3.7)  has  at  most  two  fixed  points  and  Eq.  (3.1)  allows  at 
most  two  solitary  waves.  When  Eq.  (3.7)  has  two  fixed  points,  S'(jj)  will  have  opposite  signs 
at  these  points.  From  Eq.  (3.9)  we  see  that  the  sign  of  In  a^dd  is  always  positive,  thus 
one  of  these  two  fixed  points  is  stable,  and  the  other  unstable.  We  can  also  readily  show  that, 
if  65  4-  61C5/3  <  0,  Eq.  (3.4)  has  one  unstable  fixed  point;  if  65  +  61C5/3  >  0,  it  has  none  or 
tw'o  fixed  points  (one  stable  and  the  other  one  unstable). 

2.  C5  <  0:  In  this  case,  similar  results  can  be  obtained.  Here 

[  al<W  =  y'—Z/ci  tanh"1  v.  (3.13) 

J -OO 


SM_nJh(,i+4su’+diri7)' 


(3.14) 


where  Si  (i  =  1,2,3)  are  given  in  Eq.  (3.11),  and  v  =  \/ — 16csu^/3.  Similarly,  we  can 
show  that  Eq.  (3.1)  has  at  meet  two  solitary  waves,  at  most  one  stable.  Furthermore,  if 
365  -463c5  -  167c§/3  <  0,  Eq.  (3.4)  has  one  unstable  fixed  point;  if  365  -463C5  -  167^/3  >  0. 
it  has  none  or  two  fixed  points  (one  stable  and  the  other  unstable). 
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When  Eqs.  (3.3)  and  (3.4)  allow  stable  fixed  points,  the  corresponding  solitary  wave  may  still  be 
unstable  due  to  the  non-zero  discrete  eigenvalues  moving  into  the  unstable  region  under  pertur¬ 
bations.  We  have  shown  that  when  c5  <  0.  non-zero  discrete  eigenvalues  do  not  exist,  but  when 
c-  >  0.  two  such  eigenvalues  of  opposite  sign  exist  and  are  real.  Suppose  A^  and  -A d  are  these 
two  non-zero  eigenvalues,  and  Wd  =  (vdi-  Vd2)T  and  v~d  =  °\  L'd  corresponding  eigenfunctions. 
Upon  inserting  the  perturbation  terms  of  Eq.  (3.1)  into  Eq.  (2.40c)  and  after  some  simplifications. 


we  find  that 


[m(/c3)  =  b\(m\  4-  V'2/4)  —  7  +  2m^b^  -h  3/77565, 

(3.15) 

where 

-  «&.><» 
m'  xr«(isj,  -  «&)<»  ' 

(3.16a) 

->&><» 

“  /-“.(<&  -  '&>*>  ’ 

(3.16b) 

and 

-<&)<» 

ms  jr. -  <&)■»  ' 

(3.16c) 

If  we  take  V(=  0)  and  uj  as  the  stable  fixed  points  of  their  evolution  equations  (3.3)  and  (3.4),  then 
Im(fc3)  in  Eq.  (3.15)  can  be  evaluated.  If  it  is  negative,  according  to  Eq.  (2.41),  the  solitary  wave 
in  Eq.  (3.1)  will  be  unstable. 


Lastly,  we  consider  the  instability  of  the  solitary  waves  (1.2)  in  Eq.  (3.1)  caused  by  the  continuous 
eigenvalues  of  L  under  perturbations.  In  this  case,  the  results  in  [24]  indicate  that  when  b\  <  0 
or  7  >  0.  these  continuous  eigenvalues  will  move  into  the  unstable  region  due  to  perturbations. 
Otherwise,  this  type  of  instability  is  absent. 

Now  we  summarize  the  above  results  on  the  stability  of  the  solitary  waves  (1.2)  in  the  perturbed 
cubic-quintic  nonlinear  Schrodinger  equation  (3.1).  When  c3  =  -1,  all  the  solitary  waves  are 
unstable.  When  c3  =  l  and  c5  <  0,  the  solitary  wave  is  stable  if  and  only  if  6i  >  0.  7  <  0.  V'  =  0 
and  u;  is  the  stable  fixed  point  of  Eq.  (3.4).  When  c3  =  1  and  c5  >  0.  it  is  stable  if  and  only  if 
b\  >  0.  7  <  0,  V  =  0,  uj  is  the  stable  fixed  point  of  Eq.  (3.4),  and  Im(fc3)  given  by  Eq.  (3.15)  is 
positive.  Comparison  of  these  results  with  those  by  the  adiabatic  perturbation  method  [8|  shows 
that,  when  c3  =  1  and  c5  <  0.  the  adiabatic  method  yields  the  correct  stability  conditions:  but 
when  c3  =  1  and  C5  >  0,  it  does  not.  The  reason  is  that  it  misses  the  instability  caused  by  the 
additional  mode  with  the  non-zero  discrete  eigenvalue,  Ad,  into  which  the  solitary  wave  could  emit 
energy,  as  manifested  by  Eq.  (2.41). 

As  two  examples,  we  choose  c3  =  1 ,  c5  =  1  or  - 1 ,  7  =  -0.1,  63  =  -1  and  determine  the  regions  in 
the  (61,65)  plane,  for  the  existence  of  stable  solitary  waves  by  using  the  above  results.  For  the  given 
61  and  65  values,  we  first  numerically  determine  the  stable  fixed  point  uj  from  Eqs.  (3.7)  and  (3.8) 
using  Newton’s  method.  This  alone  will  give  the  region  of  stable  solitary  waves  for  c5  =  -1.  It  is 
shown  in  Fig.  2  (II).  When  C5  =  1,  we  take  the  stable  fixed  point,  uj,  of  Eq.  (3.7),  and  numerically 
determine  the  discrete  eigenmode  A d  and  t fa  from  Eq.  (2.14),  using  the  shooting  method.  We 
then  evaluate  mk  ( k  =  1,2,3)  from  Eqs.  (3.16)  and  Im(fc3)  from  Eq.  (3.15),  with  V  being  taken 
as  zero.  The  stable  region  is  the  set  of  (61,65)  points  where  the  stable  fixed  point,  u J,  exists  and 
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Im(Jk3)  is  positive.  This  region  is  shown  in  Fig.  2  (I).  At  any  point  in  these  regions.  Eq.  (3.1)  has 
exactly  one  stable  solitary  wave.  We  observe  that  in  both  cases,  when  the  diffusive  term  ( b ■ )  in 
the  perturbation  increases,  for  the  solitary  wave  to  be  stable,  the  nonlinear  absorption  (65)  has  to 
decrease,  but  not  be  below  a  certain  lower  bound.  An  interesting  fact  is  that,  in  the  c5  =  1  case, 
stable  solitary  waves  exist  even  when  65  <  0  (see  Fig.  2  (I)).  In  this  case,  both  the  nonlinear  effects 
in  the  perturbations  are  amplifying,  but  they  are  offset  by  strong  diffusion.  Thus  a  stable  pulse  is 
still  possible.  However,  such  stable  regions  are  very  small.  Note  that  in  Fig.  2  (I),  inside  the  region 
enclosed  by  the  two  solid  curves  and  the  65  axis,  Eqs.  (3.3)  and  (3.4)  have  a  unique  stable  fixed 
point.  This  is  the  region  captured  by  the  adiabatic  perturbation  method  [8|.  But  the  solitary  wave 
is  still  unstable  in  most  of  this  region  due  to  instability  of  eigenvalue  Ad  bifurcations. 


4  Discussion 

In  this  paper,  we  studied  the  stability  and  evolution  of  the  solitary  waves  in  perturbed  generalized 
nonlinear  Schrodinger  equation  (1.1),  and  the  perturbed  cubic -quintic  nonlinear  Schrodinger  equa¬ 
tion  of  Ginsburg-Landau  type  (3.1)  in  particular.  We  found  that  the  solitary  waves  in  Eq.  (1.1)  are 
subject  to  three  types  of  instability  which  are  associated  with  the  bifurcations  of  the  zero,  non-zero 
(discrete),  and  continuous  eigenvalues  of  the  linear  operator  L  in  the  presence  of  perturbations. 
Our  stability  conditions  for  the  solitary  waves  are  both  necessary  and  sufficient.  When  specializing 
to  the  perturbed  cubic-quintic  nonlinear  Schrodinger  equation  of  Ginzburg-Landau  type,  we  proved 
that  for  any  set  of  parameters,  Eq.  (3.1)  has  at  most  one  stable  solitary  wave.  We  also  specified  the 
parameter  regions  of  stable  solitary  waves  and  graphed  them  for  two  particular  examples.  When 
compared  to  the  results  in  [8|  by  the  adiabatic  perturbation  method,  we  found  that  the  adiabatic 
method  missed  certain  types  of  instability,  especially  the  one  associated  with  the  bifurcation  of  the 
non-zero  discrete  eigenvalues  of  the  operator  L. 

The  method  we  employed  in  this  work  is  based  on  the  completeness  of  the  bounded  eigenstates  of 
the  operator  L  and  a  standard  multiple-scale  perturbation  technique.  The  key  in  this  analysis  is  the 
completeness  of  L's  bounded  eigenstates  in  L2  space.  It  allowed  us  to  solve  the  relevant  linearized 
equations  at  various  orders  and  detect  secularities  in  the  linear  solutions,  which  then  set  the  stage 
for  the  multiple-scale  perturbation  method  to  come  into  play.  For  general  perturbed  nonlinear 
wave  systems,  if  this  completeness  of  the  bounded  eigenstates  of  the  associated  linear  operator  can 
be  established,  then  the  analysis  in  this  paper  can  be  adapted  to  those  systems  as  well  and  a  full 
account  of  the  stability  and  evolution  of  permanent  waves  in  the  presence  of  perturbations  can 
be  provided.  The  completeness  of  the  bounded  eigenstates  of  a  linear  operator  has  been  studied 
extensively  in  the  literature  (see  (20,  25,  26]  for  example).  It  has  been  well  established  for  self- 
adjoint  operators.  For  generic  non-self-adjoint  operators,  as  discussed  in  the  appendix,  we  can  prove 
the  completeness  using  the  direct  scattering  technique  similar  to  that  in  [20] .  For  general  operators, 
corresponding  to  discrete  eigenvalues,  generalized  eigenstates  as  well  as  the  regular  eigenstates  may 
exist.  But  if  the  set  of  the  discrete  eigenvalues  is  finite,  we  can  still  show  that  the  eigenstates  and 
the  generalized  eigenstates  of  the  linear  operator  form  a  complete  set.  The  details  will  be  discussed 
elsewhere.  In  this  light,  our  recipe  for  the  study  of  stability  and  evolution  of  permanent  waves  in 
perturbed  nonlinear  systems  can  be  widely  applied. 
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Appendix 

[n  this  appendix,  we  study  the  spectrum  structure  of  the  operator  L  given  by  Eq.  (2.12),  and  estab¬ 
lish  the  completeness  of  its  bounded  eigenstates  in  L 2  space.  For  the  exactly  integrable  nonlinear 
Schrodinger  equation,  the  eigenstates  of  L  are  related  to  the  squared  Zakharov-Shabat  eigenstates. 
Thus  the  completeness  of  L’s  bounded  eigenfunctions  can  be  established  by  the  inverse  scattering 
technique  [15,  27].  But  for  the  generalized  nonlinear  Schrodinger  equation,  that  connection  breaks 
down.  In  this  case,  we  will  use  the  direct  scattering  method  as  developed  in  [20.  21]  to  accomplish 
this  task.  For  convenience,  we  will  replace  9  by  x.  We  first  consider  the  general  potentials  q{x)  and 
r(x)  which  vanish  at  infinity,  then  specialize  to  the  present  case  where  q  and  r  are  given  by  Eq. 
(2.8). 

The  eigenvalue  problem 

L 

can  be  written  out  as 

uxx  -  (u  +  \)u  = -r(x)u  -  q(x)v,  (A2a) 

vxx  —  (u>  —  \)v  =  -r(x) v  -  q(x)u.  (A2b) 

To  avoid  dealing  with  the  branch  cuts  at  A  =  ±u>,  we  make  the  following  parameter  transform 

A  =  u,(C2+C2)/2-  (A3) 

Then  Eq.  (A2)  becomes 

/  0  1  0  0 

_  62  —  r  0  —q  0 

Yx  ~  0  0  0  1 

^  —q  0  t/2  —  r  0 

where  Y  =  (u,ux,v,vx)T ,  and 

6  =  ^7/ 2(C  +  Cl),  r?  =  (C  -  C-1)-  (A5) 

The  rest  of  the  analysis  is  analogous  to  that  for  a  n-th  order  scalar  scattering  problem  considered 
in  [20].  Thus  the  results  will  only  be  sketched  here  with  the  proofs  omitted.  We  define  the  singular 
set  E  as 

£  =  {£  :  the  real  parts  of  any  two  numbers  of  6 ,  — <5,  rj  and  —rj  are  equal}.  (A6) 
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It  is  easy  to  see  that  E  is  the  set  of  all  the  rays  originating  from  s‘  =  0  with  angles  being  the 
multiples  of  tt/4.  We  number  these  rays  cyclically  as  Eo,  El —  and  the  sectors  C\  E  as  P.\.  Q2 — 
(shown  in  Fig.  3).  On  E0.  Re(r?)=0;  on  Eb  Re(<5)  =R e(— r?)  (r  0):  on  E2.  Re(6)=0:  in  sector  Q;. 
Re(0)>Re(-f?)>Re(^)>Re(-<5);  in  fl2.  R e(-r?)>  Re(<5)>Re(-<5)>Re(r?):  etc.  In  each  sector,  we 
define  two  fundamental  matrices  and  of  Eq.  (A4)  according  to  the  ordering  of  6.  -6.11  and 
— 1 ]  in  that  sector.  For  instance,  in  fl\.  we  define 


<&“(x.  £)  =  m~(x,  OelJ ■ 


where  J  =  diag(6.  —t),  rj ,  -<5), 

m~{x,  C,) 


(  1  0  0  1  \ 

<5  0  0  — <5 

0  110 
V  0  -17  V  0  ) 


X  — >  dbOC, 


1  AT) 


(A8) 


and  m~  are  bounded  as  |x|  — >  oc.  In  other  sectors.  <t>~  can  be  similarly  defined.  It  is  easy  to  see 
that  <&-  so  defined  are  unique.  They  exist  for  all  £  €  C  \  E.  apart  from  a  discrete  set  Z  which  are 
all  the  zeros  of  A*  ( k  =  1,2,  3)  to  be  defined  below.  At  the  points  Cfc  €  Z,  have  pole  singularity. 


Next  we  define  the  functions 

Ai  —  mi  A  m2  A  m3  A  /( — 4677) ,  (A9a) 

A2  =  771 3*  A  7712  ^  m3  ^  m4  /(“ (A9b) 

A3  =  mj"  A  TTij  A  7713  A  7714  /(— 4<5p),  (A9c) 

where  “A”  represents  the  wedge  products  of  vectors.  It  can  be  shown  that  Ak  ( k  =  1.2.3)  are 

independent  of  x,  analytic  in  each  sector  0*,  and  Ak  —  1  as  |Cl  —  00.  Furthermore.  A2  is  analytic 
across  the  boundary  E 1 .  The  discrete  eigenvalues  A  of  the  operator  L  correspond  to  the  zeros  of 
A2  through  relation  (A3).  The  continuous  eigenvalues  of  L  correspond  to  the  two  rays  Eo  and 
E2.  We  choose  a  path  P  in  the  £  plane  as  shown  in  Fig.  3,  with  its  direction  counter-clockwise. 
This  path  starts  at  £  =  0+  +  ooi,  moves  down  vertically  to  £  =  i,  half-circles  around  it.  and 
moves  downward  again  until  it  reaches  £  =  0.  Then  it  quarter-circles  around  £  =  0  and  moves 
horizontally  along  the  upper  side  of  the  positive  Re(£)  axis,  until  it  arrives  at  £  =  1.  Then  it 
half-circles  around  £  -  1,  keeps  on  moving  horizontally,  and  eventually  ends  at  £  =  oc  +  0~j.  In 
the  A  plane,  this  path  corresponds  to  one  which  encloses  the  entire  A  plane  except  the  continuous 
spectrum  {A  :  A  >  u  or  A  <  -u>}.  Thus  the  winding  number  of  A2  along  P: 

N  =  d-r  [  TTTTrdC  =  ^"{arg{A2(oo)  -  arg{A2(co  +  i)}  (A10) 

2irt  Jp  A2K) 

gives  the  total  number  of  the  discrete  eigenvalues  of  L  (multiplicity  of  non-simple  eigenvalues 
counted). 

The  completeness  of  the  bounded  eigenstates  of  L  in  L<i  space  can  be  established  by  constructing 
the  Green’s  function  to  the  equation 

(L  -  A)G(x,  y,Q  =  «5y(x)diag(l,  -1),  (All) 
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and  proving  that  G  can  be  expressed  as  a  linear  combination  of  L's  bounded  eigenstates.  We  call 
the  operator  L  generic  if  (1)  A k  (fc  =  1-2.3)  have  no  common  zeros  and  no  multiple  zeros.  (2) 
they  have  no  zeros  on  A  and  (3)  the  set  of  their  zeros  is  finite.  For  a  generic  self-adjoint  operator, 
the  completeness  of  its  bounded  eigenstates  in  L-2  space  was  proved  in  [20]  using  this  approach.  In 
our  case.  L  is  not  self-adjoint.  But  if  it  is  generic,  slight  modification  to  the  analysis  in  [20]  can 
be  made  to  establish  the  completeness  relation  as  well.  If  L  is  non-generic  (as  is  the  case  when 
q(x)  and  r(x)  are  given  by  Eq.  (2.8)).  a  discrete  eigenvalue  may  be  multi-fold,  and  its  algebraic 
multiplicity  may  be  larger  than  its  geometric  multiplicity.  In  such  cases,  the  generalized  eigenstates 
of  the  discrete  eigenvalues  are  also  needed  for  closure  (see  [15]).  But  as  long  as  the  number  of  L  s 
discrete  eigenvalues  is  finite,  similar  techniques  can  be  used  to  show  that  the  bounded  eigenstates 
of  L  (including  the  generalized  discrete  eigenstates)  also  form  a  complete  set.  As  a  practical  guide, 
if  a  discrete  eigenvalue  is  a  k- th  fold  root  of  A2,  then  correspondingly  k  regular  or  generalized 
eigenstates  should  be  included. 
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A 


Figure  1:  The  trajectory  of  A2(0  as  C  moves  along  the  path  P  for  C3  =  1  and  C5  <  0.  The  path  P 
is  specified  in  the  appendix  and  shown  in  Fig.  3. 
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(I) 


(II) 


Figure  2:  Stability  regions  (shaded)  of  the  solitary  waves  in  Eq.  (3.1)  for  C3  —  1,7  =  —0.1  and 
63  =  -1.  (I)  C5  =  1;  (II)  C5  =  -1. 
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Figure  3:  The  complex  £-plane. 
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